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Abstract. A simple and natural explanation for the dynamics and mor- 
phology of the Local Group Dwarf Spheroidal galaxies, Draco (Dra) and 
Ursa Minor (UMi), is that they are weakly unbound stellar systems with 
J> ■ no significant dark matter component. A gentle, but persistent, Milky 

r**"- \ Way (MW) tide has left them in their current kinematic and morphologi- 

cal state. This short paper reviews the parametric tidal interaction model 
\ which accounts for their behavior and discusses new statistical evidence 

■ from the observed stellar distribution in Dra which implies that its total 

ff"} \ mass is not dominated by collisionless dark matter (DM). 
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The Milky Way may be our best testing ground for models of cosmologically 
important dark matter masses. The dumpiness of the MW potential and the 
mass and number density of dwarf (M< 1O 8 M ) galaxies is an important con- 
straint on "standard" Cold Dark Matter (CDM) mass components (Mayer et al. 
2002). The paucity of dwarf galaxies in the Local Group is something of a crisis 
for CDM (Klypin et al. 1999) and adds to the urgency of understanding what 
the masses of the MW dS really are. 

Controversy over the dynamical masses of the MW dS centers around 
whether or not they are in equilibrium. At one extreme the equilibrium models 
account for elliptical morphology and dynamics with ad hoc multiparametric 
dark matter mass distributions (cf. Kleyna et al. 2002). In contrast a variety 
of non-equilibrium dS models account for their properties through MW tidal 
interactions (Kuhn and Miller 1989; Kroupa 1997). Several years ago Kuhn and 
Miller (1989) argued that if only the dS luminous mass component is present (as 
is the case in globular clusters) then their dynamical timescales are comparable 
to their MW orbital times. This has important consequences since then even a 
relatively weak MW tidal force can profoundly affect the dS through its orbital 
time dependence. Recently we developed analytic and numeric models for a 
parametric resonance that describes how a variable tidal interaction inflates the 
velocity dispersion of a dS stellar system (Fleck and Kuhn 2003 - FK). 
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In a parallel empirical effort we have been looking for tidal tails as evidence 
of dS non-equilibrium. The Sagitarrius dS tail is the most obvious example of a 
dS-MW tidal interraction which, incidently, was discovered (Ibata et al. 1994) 
after the early non-equilibrium tidal predictions. Other examples of tails beyond 
a dS tidal radius and with a surface density far above what is expected from a 
massive equilibrium system have been discovered in Carina (Kuhn et al. 1996), 
Ursa Minor (Smith et al. 1997), and Sculptor (Walcher et al. 2003). 

Here we also briefly describe a new test of the simplest equilibrium dark 
matter hypothesis in dS using relatively deep color-magnitude photometry from 
our CFHT wide-field survey of Ursa Minor and Draco. By color-selecting and 
then counting stars into the dS main sequence the statistical distribution of the 
stellar number density becomes a sensitive probe of the overall dS gravitational 
potential. 

2. Parametric Tidal Oscillations 

The cumulative effect of a weak MW tide acting on a dS has been described by 
FK in terms of a modified two-dimensional oscillator. This equation is analogous 
to a one-dimensional Mathieu equation (cf. Bender and Orszag 1978) which 
describes an inverted pendulum or the amplitude of a "pumped" swing oscillator. 



The dS-centered local cartesian coordinates describe the position of a star 
and ip is the orbital phase angle of the dS in its MW orbit. Here e = oJI/ooq 
is the strength of the tide expressed in terms of the ratio of the MW circular 
frequency and the harmonic frequency of the self-gravitating dS. Solutions to 
this equation are known as parametric oscillations and have several interesting 
properties. 

Unlike driven simple harmonic oscillations, where the growth rate of the 
oscillation depends strongly on the intrinsic mode damping, parametric modes 
can exhibit secular growth even when the intrinsic harmonic oscillation frequency 
is significantly different from the driving frequency. If the orbital phase of the dS 
is ip = tot then FK showed that 2-d parametric modes are unstable when loq{\ — 
e/4) < lo < oJo(l+e/4). Non-harmonic dS potentials exhibit qualitatively similar 
behavior. FK found that a large fraction of possible dS orbits are subject to 
parametric excitation. They also showed that parametric oscillations generally 
cause an inclined dS bar with an inclination angle that depends on the dS-MW 
orbit and potential, and with a growth rate which is a simple function of e. 

From numeric and analytic solutions to the tide equations FK showed that 
centrally concentrated dS stellar distributions, but with inflated velocity disper- 
sions, could be generated after several MW orbits. In general the tidally induced 
velocity dispersion of any dS system should increase with s. This dependence is 
consistent with the observed MW dS. Figure 1 shows how the "apparent M/L" 
varies with the empirical mode growth estimate (A = yfe) based on a standard 
logarithmic MW potential and the luminous mass of each of the dS. The large 



x + ojq 2 (1 — e cos 2ip)x = eu!o 2 ysin2ip 
y + ooq 2 (1 + e cos 2tp)y = eu;o 2 xsin2<£ 
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(2) 



APS Conf. Ser. Style 



3 



2.5 7 
2.0 - 

27 1.5- 
cn 

23 1.0 : 

0.5 - 
0.0 2 

-1.2 - -- - -0.4 

Log (A) 

Figure 1. Variation of dS apparent M/L with parametric mode 
growth estimate. 

apparent M/L ratio of each of the local dS is consistent with this model of 
parametric tidal excitation, leaving no kinematic evidence for a significant dark 
matter binding mass. 

3. Stellar Number Density Statistics 

There exist equilibrium mass calculations for Dra which claim a M/L in excess 
of 500 solar units (Kleyna et al. 2002). With this assumption the Dra potential 
must be dominated by its dark matter. If we assume the DM is collisionless 
and relaxed, and that it and the baryonic mass of the dS formed at the same 
epoch, then the dS stars are, effectively, massless test particles within a smooth 
DM-dominated potential. Such stellar "particles" must be uncorrelated, with a 
null two-particle correlation function (cf. Binney and Tremaine 1987). It follows 
that the stellar particle distribution function, f(f,v,t), can be computed from 
the collionless Boltzmann equation and the net dS gravitational potential (cf. 
Binney and Tremaine 1987). 

For a smooth gravitational potential we also expect f(f) to be smooth so 
that the projected star counts in a small spatial region (for example as derived 
from one pixel of a two-dimensional star count map) can yield a reasonable 
estimate of the projected number density, n(x). Here x is a 2-d coordinate in the 
plane of the sky and a line-of-sight and velocity integration of f(f, v, t) to obtain 
n(x) is implicit. On scales small compared to the variation of the potential we 
also expect n(x) to be smooth. In fact, since stars are uncorrelated, our estimate 
of n(x) from the star counts per pixel should be Poisson distributed - like the 
photon count in a 2-dimensional image. Any spatial "dumpiness" in the number 
of stars per pixel which is in excess of Poisson distribution variations violates 
our assumption that Dra has a smooth DM dominated potential. 

Using deep V and I color-magnitude data obtained from the CFHT 12K 
CCD mosaic camera over a 2 x 1° region centered on Dra (Kocevski et al. 
2003) we estimate the circularly symmetric projected number density n{R) when 
R = \\x\\ (Figure 2). The effect of non-dS background stars was minimized by 
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Figure 2. The number density of Dra stars in elliptical bins versus 
projected distance from the Dra center. 

counting stars in a color-magnitude region that included the Dra giant branch, 
horizontal branch and main sequence down to V=25. Radial bins were scaled 
along the y axis to account for the ellipticity of Dra. The radial binwidth also 
increased with R in order to maintain a constant count per bin, and constant 
statistical significance, since the mean density n(R) decreased with R. 

Figure 2 shows n(R) in approximately 10 3 radial bins from a region slightly 
larger than 1° in radius centered around Dra. The fluctuations in n{R) in groups 
of 5 bins were used to compute the normalized (reduced) chisquared. At radial 
bin Ri we obtain (with the assumption of Poisson statistics) 

I 4 

5 i=o 

Here <7 2 (i?j) is the sample variance computed at each radius bin from it and 
4 nearby bins. We find that x 2 derived from the outer bins (with R > 20') is 
consistent with the expected reduced \ 2 distribution. The inner Dra bins have 
a variance which is too large. Figure 3 plots the distribution of x 2 measured this 
way for the inner Dra region. Using a Kolmogorov-Smirnoff test we find that 
the inner region is not consistent with the expected x\ distribution at better 
than a 99% confidence level (while the outer region is). Thus, the stellar density 
in the inner region of Dra is "clumpier" and more variable with spatial position 
than it should be if the Dra mass were dominated by a smooth DM component. 
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Figure 3. The actual (solid) and expected (dotted) x| for the stellar 
number density variance in the inner 20' of Dra is plotted here. The 
preponderance of larger-than-expected empirical %| values implies that 
the potential of Dra is not dominated by a smooth DM mass distribu- 
tion. 
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